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The present invention relates to a system for analysing ECG curvature where at least 
one among a number of dif&ient parameters is isolated, which system has a input 
means connected to an ECG source, where the difif^^t parameters of a received ECG 
curvature are indicated and/or isolated and for indicating possible symptoms which 
S relates to or are indications of certain deceases, where said deceases are known to m-* 
fluence the ECG curvature. 



The present invention further relates to a method for analysing ECG curvature, which 
curvature contains a number of parameters. 

The aim of the invention is to achieve a system and a method for objective^ &st and 
effective indication of a number of symptoms derivable fiom an ECG curve which 
may be indicative of one or more diseases. 



1 5 This can be achieved with the system previoxisly described if a first number of selected 
parameters, are combined in at least a first mathemadcal analysis, where the result of 
the analysis can be represented as a point in a coordinate system comprising at least 
two axes where the system can compare the actual placement in the coordinate syst^ 
with a number of reference parameters stored in the system for indicating diseases 

20 having influence on Ae ECG curvature. 



Hereby,- it is achieved that any symptom of a disease having an indicajSon (inflomice) ~ 

in the ECG curvature can be detected in an objecnve, automated and veiy fast way. 
The system might be used under field conditions such as in ambulances or in other 

25 situations where a fast indication of heart diseases is needed in order to help the pa- 
tient in a correct way as early as possible. The analysis that takes place in an ambu* 
lance on its way to the hospital, can by transmitting the residts to the hospital, allow 
the doctor at the hospital to give feedback to the personnel in the ambulance so that 
the correct treatment of die patient may start At the same time, the hospital can pie-* 

30 pare the correct activity for the incoming palieat. Hie system could be very important 
. for ECG analyses for all non-specialists in the field if ttiey have to analyse a ECG 
curvature. 
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Hie scope of the inventioa can also be fulfilled with a method for analjrsing the ECO 
curvature if the method incorporates the steps of: 



- leceiving BCG curvature from a source^ 

5 - indicating a number of different parameters contained in the received ECG 

curvature, 

^ storing the parometeis in storage means^ 

" selecting disease specific parameters in the storage mean$ 

^ combining selected parameters in mathematicaL analysing meaiis 

10 - representing the result of the ma&ematical analysis as ft poin* i^ 

system, comprising at least two axes» 

- comparing the actual placement in the coordinate system with a manber of ref- 
erence parameters stored in a memory, 

- indicating possible diseases having the determined influence on the ECG cur- 
15 vature. 

In this as aheady described, a v^y effective analysis of the ECG curvature is 
achieved. 



Below the parameters are illustrated on a typical ECO curvature in arder to illustra 
thejiiSerent curve sections isolated by tbe analysis, and refeired to by parametw^ 
intovals aiu! complexes. 



Pmo* 7«it:i.q/l7/!?003 11:06 Empf.nr*:897 P.011 



NO- 310 



3 





JO 



15 



20 



The heart generates an electrical signal ECG. The waves at the ECO-signal P OJLS T 
andUareduetodepolarisationaijdrepolflrisaiionofthehean. ' 



T^QRSKTompIex contains the Q^R and S wave and goes &om the Q^mset to S- 
The QT interval starts at the beginning of QRS (Q-onset) and ends at the T- 



■>«ave on- 



TTiB RR interval goes from one R-peak to the following. 
The PR interval is from the beginning of P to the beginning of QRS (Q-onsei). 
QRS duration is the width the QRS complex. 
The ST-segment starts at the QRS-complex offset and ends at the T- 



wave onset 
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The analysing process is repeated in the system for further selected parametexs in order 
to indicate fiirtber possible diseases or symptoms. Hmby^ it is achieved that the sys- 
tem or Che method can be repeated several times with difierent combinaiiDns of pa- 
'\ ~ rameters. In order to analyse for a high number of possible diseases^ each time a result 

5 of an analysis is adiievedj, the parameters are compared with stored parameters where 
eveiy stored paramet^ is an indication of a symptom of a known disease. Witti the 
system, a deviation of parameters from tbe stored data indicating symptoms of an ex- 
act disease may also be inteipreted for further reference. Parameters could be con- 
trolled, and if ^e storage means has an enormously high number of reference paiame- 
10 . ters relating to vanons symptoms^ testing fi>r all symptoms and correlated diseases 
having influence on ^ ECG curvature may be peiformed. 

The system or method divides the paramet^a into at least two main groups^ which 
groups contain paEameters of symmetiy, flatness and duiatioii relating to the actual 
15 ECO curvature. In this way» it is achieved that the paranieters are grouped in Oie sys- 
tem, and in every gcoup^ they can be further subdivided into a specific nmnber of pos- 
sible parameters. Keeping the number of parameters relatively small, Ae analysis takes 
place in a faster way. 

20 The group of symmetry might contain at least the following parameters: 

5 1 Skewness evaluated from Tstart to Tend, 

* 

52 Skewness evaluated fix)m Tstart to Tend with Ttop as mean, 

53 Skewness evaluated in a symmetric interval, 10 % of the Tstart-^Tend interval 
suiTounding Ttop with Ttop as mean 

25 S4 Skewness evaluated in a symmetric interval, 20 % of the Tstait-Tend interval 

sunounding Ttop wit Ttop as mean, 

S5 RatzD of the time int^al firom Tstart to Ttop and the time interval from Ttop to 
Tend. 



S6 Ratio of the average slope j&om Tstart to Ttop and from Ttop to Tend. 



30 
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The group of flatness might contain at least the foUowing paiameteis: 



. —El — Kuttosds evaluated fixnn Tatait to-Toadr 



Fl nomielized by thie absolute Rtop-Qnadir value, 

5 F3 Kurtosis evaluated from Tstait to Tend with Ttop as mean, 
F4 F3 nonnalized by absolute Rtop-Qnadir value, 

F5 Kurtosis evaluated in a symmetric interval. 10 % of the Tstwt-Tend interval sur- 
rounding Ttop with Ttop as mean. 

F6 F5 noraaalized by absolute Rtop-Qnadir valuer 

10 F7 Kurtosis evaluated m a jg^nmieiric internal, 20% of the Tfl^^^ 

rounding Ttop with Ttop as mean, 

F8 Kurtosis nonnalized by the value of Rtcp with Ttop as mean, 

F9 Ratio of the total area under the T-wave from Tstart to Ttop and the coire^ond- 
xng time inteirvdl, 

15 FIO F9 nonnalized by absolute Rtop-Qnadir value, 

FlI Ratio of the total axea under the T-wavB from Ttop to Tend and the conespond- 
ing time interval, 

Fl 2 Fl 1 nonnalized by absolute Rtop-Qnadir value 

F13 .Mo of the total area under-theT-wavefrom-Tstartto-Tend and the cwrespoad- 

^ ing time interval, 

F14 F13 nonnalized by absolute Rtop-Qnadir value. 
F15 Ratio of the height of Rtop and the widflj of the Tsian-Tend interval. 



Hie group of duration might contain at least the following paiameteis: 



QTc The Q-T imerval nonnalized by the squareroot of the R-R interval accorfiiig to 
Bflzett's fonnula , 

D2 Time interval from Tstart to Teai, 
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D3 Time interval fiom Tstart to Ttop, 
D4 Time interval frova Ttop to Tend. 

The groups of parameters could contain further paiameteis and die group may contain 
5 a numbCT of sub groups. 

Combining parameters, the system aod/or method may form the different groups. 
When combining parameters from different groups, a much better result is achieved 
than when only using parameters from the same group* The parameter? can be an ele- 
10 vation of the curve; they can be the moxpholcgy of the curve; or they could be lime- 
deviations as an example of possible parameters. When combining paiameteis, a pre* 
dse analysis can take place because a specific combination of parameteis can indicate 
a specific disease and it is possible effectively to select between ECG*signals that look 
very much alike, but which indicate different diseases* 

IS 

The system and/or method can analyse the QT interval of the ECG curvature fi>r indi- 
cating Long QT syndrome. This w^yj the Long QT syndrome can be indicated in an 
objective and effective manner which might occur v^th children right after they are 
bom in order to start a treatment of Long QT syndrome as early as possible. The 
20 method can differentiate between different genotypes of the Long QT Syndrome, 
vAdch is important for the treatment Hereby can be achieved that the correct medical 
treatments can be started. 

ITie system and/or method can analyse fi>r ST-elevation myocardial infarction by anar^ 
25 lysing ax least the following parameter?: ST elevation, ST morphology, T wave: mor- 
pholep and Q wave moiphology* In this way, a very effective indication for ST- 
elevation myocardial infirction is achieved at the correct activity can be stated as early 
as possible. 

30 The system and/or method caa analyse for Non ST-elevation myocardial in&rction by 
analysing at least the fbllowii^ parameters: ST depression, T wave moiphology and Q 
wave morphology. Non ST-elevation myocardial infiarction can also be detected in a 
highly effective way and correct treatment can be stated* 
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The system and/or method caa analyse for Cardiomyopathia by analysing at lease ihe 
followmg paiameters: P vmve moiphology, QRS duration, S Wave morphology, T 
wave morphology. Cardi^omyopathia can in this way be effectively detected. 

5 

The system and/or method can analyse for Brugada Syndrome by analysing at least the 
following parameters: PR^duration, ST elevation^ ST morphology and T wave mor- 
phology. Hereby, it is achieved that Brugada Syndrome is detected in a veiy effeictive 
and fast way. 

10 

Ihe system and/or method can analyse for Right bundle branch block KBBB by ana- 
lysing at least the following parameters: QRS duration, QRS morphology, T wave 
morphology and ST elevation. In this way, an indication of Right bundle branch block 

■ 

RBBB is achieved in an effeciive way. 

15 

The system and/or method can analyse for Left bundle branch block LBBB by analys- 
ing at least the following paramelBrs: QRS duration, R wave moiphology and T wave 
morphology. An efifecdve indication of Left bundle branch block LBBB is also possi- 
ble with this method or by using the system. 

20 

The ^y3tem and/o r mediod ca a analyse for Short QT Syndromejby jiialysing.a^ 
the following parameters: Q-T dviation and T wave morphology. Short QT Syndrome 
can also be analysed effectively. 

23 The system and/or method can analyse for Hyperkalemia by analysing at least the fol- 
lowing parameters: P wave morphology, T wave moiphologyf QRS duration^ QT dura- 
tion and PR duration. Hereby. indicatioA of Hyperkalemia is achieved in a highly ef-* 
fective way. 

30 The system and/or method can analyse for Hypokalemia by aaalysiDg ar least the fol- 
loMfing parametecs: QT duration, T wave moiphology and ST depression. Hereby, it is 
achieved that Hypokalemia is indicated effectively. 
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The system and/or method can analyse for peri/myocardiiis by analysing at least (he 
follpwing parameteis: ST elevation, ST moiphology, Q wave moxphologjr and PR 
depression. Hereby, it ia achieved that peri/myocarditis is detected egeciively, 

5 The system and/or method can analyse for Right Ventricular Hypertrophy by analysing 
at least tibe following parameters: Q wave morphology, QRS duration, S wave mot- 
phology and T wave morphology. Herby. it is achieved that Right Ventricular Hyper- 
trophy is indicated efiectively. 

10 The system and/or method can analyse for Left Veatricular Hyperttophy by analysing 
at least the fbllo^ving paiametBcs: Q wave morphology, QRS duiation, S wave mor- 
phiology and T wave nuwphology. Also Left Venticular Hypertrophy can be detected 
effectively. 

1 5 The system and/or me^od can analyse for ArdvyttunogcDic Right Vwrnicular Dj^la- 
sia by analysing at least the fbllovdng patamcteis: QRS duiafioo, S wave morphology 
and T wave morphology. Hereby, it is achieved that Anhydunogenic Right Ventricu- 
lar Dysplasia is detected in an elective way. 

20 Below is described one possible method and a system to illustrate the invention. 
Abstract " 

The Long QT Syndrome is a genetic disorder characterized by abnormal cardiac repo- 
laiisatioii resulting m prolonged QT duration, syncopal episodes and increased risk of 
sudden cardiac death. Mutations in the KvLQTl- and HERG genes account for more 

25 than 90% of all LQTS patients. The^ interval duration is the only ECQ-based quan- 
tifier of LQTS used in dinical practice today. However duiation is only a gpjss esti- 
maie of lepolarisadon and does not allow pofect discrimination between KvLQTl, 
HERG and normal subjects. Studies have shown that T-v««ve morphology parameters 
are useful discriminators in LQTS, but no single parameter has proven to be sufficient. 

30 la this study we present a novel multivariate disodminatiQn method based on a combi- 
nation of T-wave symmetry-, flatness- and duration parameters. 16 subjects were in- 
cluded m the shidy - 8 normal, 5 HERG and 3 KvLQTl patients. Genotypes were 

P„o* ^5t-iq/i7/5nQ3 11:07 Enpf.nP.;:897 PJ3.17 • • 
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known for all LQTS patients^ but ooe. Standwl 12 - lead ECG's were recorded on 
each siihject An automatic ECG event detection algorithm was implemented. The 
signal was highpas5 filtered and noimalized with respect to the isoelectric level to en- 
sure a stable baseline. 4 parameters describing the duratioa of repolarisations 6 syzn- 
S metiy- and 15 flatness pai^meteis were calculated to characterize each of the T-waves. 
The mean values of lead VS and the interlead standard deviations were used as pa- 
rameter values* Stepwise discriminant analysis was performed to obtain two discrimi- 
nant iimctLons based on the five strongest discriminatoiy parameters. Tlie resulting 
discrindnant functions include Z duratton^» 2 symmetiy- and 1 flatness parameter. The 
1 0 two fiuictions classiQf all subjects correctly (p> 0,000 1 , pO.005), Further discrimlndnt 
analysis with a reduced number of parameter cate^ries implied that superior classifi- 
cation is obtained when using all diree parameter categories presented. A combination 
of parameters from the three categories symmebyy flatness and duration of repolarisa- 
tion was sufficient to correctly classify ECG recordings from the KvLQTl , HERO and 
15 normal subjects in this study. This multivariate approach naay prove to be a powerful 
clinical tool. 

1. Introduction 

Hie Long QT Syndrome (LQTS) represents a hereditary genetic disorder characterized 
by the presence of prolonged QT duration on the BCG, syncopal episodes due to po- 
lymorphic ventricular tachycardia, (torsade de pointes), and onythmogenic sudden car* 
diac death. 

* 

Mutations involving 6 different gen^ have been identified in LQTS subjects* These 
mutations result in structural and functional changes in ion-chaxmel proteins and cur- 
rents. The changes are manifest by QT prolongation and morphological gene-specific 
repolarisation patterns. The most prevalent genes affected in LQTS patients are 
KvLQTl and HERO which account for more Hm 90% of LQTS genotype patients. 
The current smdy focuses on carriers of these two genes. Although some attempts 
have been made to develop quantitative measures that link different repolarisation 
abnormalities to specific LQTS related channel-opatbities these methods have so &r 
&iled to provide a solid diagnostic yield. lo current practice the duration of the QT 
interval is the only widely accepted quantifiw of ventricular repolarisation. Yet, it has 
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been recognized that the duration of the QT interval is only a gross estimate of repo- 
larisation since T-wave morphology is also important when characterizing the QT in- 
terval. This is evidenced by the &ct that approximately lO^o of all mutafion oaniers 
have a nonnal Bazett cozrected QTc (<440md) and 40% of KvLQTl and HERG cairi- 
5 CIS shov; QTc values between 410-470 ms that overli^ with non-cani^s. Conversely 
only 2% of all carriers present with a normal ST-T pattern and a normal QT intervaL 
Morphological aberrations thus canry major implications for the identification of ab- 
normal repolarisation and have been included as diagnostic criteria equivalent to that 
of a positive &mily history for LQTS. 

10 

Studies hove shown diat affected KvLQTl patients generally show broad based T- 
waves with a normal to relatively high amplitude and often widiout a distinct T - wave 
onset For individuals with mutations involving die HERO gene flie aforementioned 
studies have g^erally found low amplitude T-waves vvith bifid T-waves in 60% or 
15 more of the carriers. 

Cardiologists already include a qualitative assessment of T-wave morphology firom the 
ECG in order to obtain information that augments the clinically established QT inter- 
val m^urwent and facilitates disprimination between LQTS genotypes. However 
20 qualitative desmption of repolarisation morphology may be biased due to intra-* and 
interpersonal variability thus indicating the need for a standardized quantitative meas- 
ure of this parameter. 

In the following is presented a novel multivariate categorization method that allows 
25 discrimination between KvLQTl, HERG and normal individuals based on Twave 
morphology recorded ftom 124ead ECG's. Halhnark morphological features of T- 
waves reported in the literature for these three groups served as inspiration for select- 
ing three primary T-wave characteristics to be assessed. These charBcteristics are 
symmetry, flatness and duration. 
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2. Methods 

2.1 Subjects 

The study uicluded ECG iieoordinga from 8 female aad 8 male sut^ects. Hie subjects 
were divided iniD four groups; 3 KvLQTl (aged 20-48, 2 females), 5 HERG (ag^ 13- 
5 76, 2 femalesX 8 nonnal (aged 23-31, 4 females). Genotypes were known for all 
KvLQTl and HERG subjects with a single exception: 1 patient was categorized as a 
KvLQTl subject by anamnesis and ECG-analysis. In the nonnal group there were no 
reports of prior cardiac diseases or LQTS fiunily pieced^t 

2.2 Data collection 

10 Data acquisition was carried out with the subjects resting in supine position. The 
equipment used for data acquisition was a portable digital ECG recording system, 
'HUardio Perfect Resting ECO system" manuiactured by Cardiocontrol. Recording was 
divided usio three sessions. Data was collected ftom 8 leads (I-m, V2-V6) with a 
sampling rate of 1200Hz. Signal recording length was 75 s. in the first session and 150 

15 s. in the last two sessions. 



Following data acquisition, SCP files generated by the Cardio perfect software were 
exported from a MSDE/SQL7 server and subsequently converted to .MAT files using 
SCP-Batch Converter. 

2.3 Algorithm for detection of events in the ECG 

To facilitate evaluation of the repolaiisation process and the QT interval, several 
events in the ECG were detected (Qstait, Rtop, Tstart» Ttop and Tend). An algorithm 
for detecting these events was implemented in Matiab 6.0. 

25 

The method is based on prior work published by Laguna et al. and uses adaptive 
thresholding techniques applied to a digitally filtered and di£^entiated signal. A mi- 
nor extension to the algorithm w£ks incorporated to enable the detection of TstacL 
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Tstart was detected with a tedmique equivalent to the lechiuqTie fox detecting Tend, 
Ficure 1 shows an ex;ample of the result of the event detection algoritkiu 



Figure 1. Importani events that are used to describe r^olarisation are marked by dots 
5 by the event detection algorithm- The algprithm is able to detect the evCTts on all 8 
recorded leads. 

2> 4 Preliminary signal processing 

Evaluation of the QT interval and the repdlarisation process was done on the basis of 
an ECG signal with stabilized baseline. This was achieved through prelirainaiy sign&l 

10 processmg- The ''taV ECG was filtered by a minimum order equuipple high pass 
filter with a cut-ofF fisquency of 2 Hz, 60 dB damping in 0 Hz and IdB tipple in the 
passband CN=1276X After filtering, the signal had an ahnosi stable baseline, &i order 
to improve stability, isoelectric lines m the signal were estimated &om one P-Q inter- 
val (Qstart minus 20 ms) to the following P-Q interval (Qstait minus 20 ms). The sig- 

15 nal was then noimalized by subtracting the line value fiom the corresponding signal 
values. This process is shown in figure 2. 

2, S T*wavc morphology parameters 

In order to chaiacteiize the T-wave morphology, a number of parameters were se- 
lected. The parameters were chosen to cover each of the three categories: Twave 
20 symmeciy, T-wave flatness and duration. The parameters are listed and described in 
table 1. 

PaianieterB S1-S4 and FI-F8 is based on tiie calculation of modified skewness and 
kuitosis measures. Inspired by the suminaiy measures of probabUity distributions used 
25 in the field of statistics the T-wavea were modelled as probability mass distributions 
(figure 3) and assigned a ceaitte (mean), widflt (standard deviation^ an asynmietry 
measiwe and a convexity measure. Asymmetry and convexity calculations we» then 
carried out based on Ae modified skewness and kurtosis measures (3'" and 4* order 
moments) as follows: 



The total area under the signal, mO, was calculated: 
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^Tbe signal was nonnalized by the value of the area; mOr"**" ' 

5 

Normalization facilitaied the calculation of the moinem functions, since w[n] shares a 
fundamenial property with fhe probability mass functioa: A totd area of 1. 
The 1^ order moment, ml, was calculated, ml is the mean of the sigoa]: 

The 2"^ order mament. ni2» was calculated, m2 is the standard deviation of the signal; 

15 

Figure 2 Isoelectric lines (dashed lines) in the signal are calculated fiom one P-Q in- 
lerval to the following P-Q interval (Qstart - 20 ms)* The Una values are subuacted 
firoqi the corresponding ECG signal values giving the distances v(n). The resuli of this 
, , . procedure is shown as an area plot with basis on the z^o-line, — — — • 

20 Figure 3. a) Exanqple probability mass distributiou used whOT calculating standard 
skewness and kurtosis measures, b) Modified ixequency distribution used in this study 
for calculating the modified skewness and kuitcsis measures. Signal values v(n) are 
shown in figure 2, 

4/8 

25 PARAMETER DESCRIPTION 

Symmetry 

S 1 Skevmess evaluated from Tstan to Tend. 
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52 Skewness evaluated from Tstart to Tend with Tlop as mean. 

53 Skewness evaluated in a symmetric inietval, 10 % of The Tstait-Tend interval 
surrounding Ttop with Ttop as mean, _ ■ ' 

54 Skewness evaluated in a symmetric interval 20 % of the Tstait-Tend interval 
5 sunounding Ttop wit Ttop as mean. 

55 Ratio of the time iniermlftom Tstart to Ttop and die lime interval from Ttop to 

Tend. 

56 Ratio of the avraage slope bom Tstait to Ttop and from Ttop to Tend. 
Flatncas 

10 Fl Kurtosis evaluated firom Tstan to Tend. 

F2 Fl xioraialized by the absolute Rtop-Qnadir value, 

F3 Kunosis evaluated from Tstart to Tend with Ttop as mean. 

F4 F3 nonnalized by absolute Rtop-Qnadir value. 

F5 Kurtosis evaluated in a symmetric interval, 10 % of the Tsiart-Tend interval sur- 
15 rounding Ttop with Ttop as mean. 

F6 F5 nomialized by absolute Rtop-Qnadir value. 

F7 Kurtosis evaluated in a symmetric interval, 20 % of flie Tstart-Trad interval sur- 
rounding Ttop wi& Ttop as mean. 
F8 . Kurtosis normalizsed by the value of Rtop with Ttop as 



20 F9 Ratio of the total area under the T-wave from Tstart to Ttop and the conespond- 

ing time interval. 
FIO F9 normalized by absolute Kiqp-Qnadij: value. 

Fl I Ratio of the total area under the T-wave from Ttop to Tend and the coirespond-^ 
ing time interval. 

25 F12 Fl 1 normalized by absolute Rtop-Qnadir valae. 

F13 Ratio of the total area under the T-wave from Tstart to Tend and the correspond- 
ing time interval. 
F14 F13 normalized by absolute Rtop-Qnadir value. 
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F15 Ratio of the height of Rtop and the width of the Tstait-Tend interval 
Dnration 

QTc Ihe Q-T interval nonnalizcd by Uie square root of the R-R inteival according 
BazBtt's formula. 

D2 Time interval from Tstait to Tend. 
D3 Time interval from Tstait to Ttop. 
134 Time interval from Ttop to Tend. 



to 



■wave 



The table above shows a Complete list of the parameters used to characterize T-^ 
morphology. Parameters belong to one of iJttee categories: symmetiy, flatness and 
duration. 



The 3 order moment, m3, was calculated. m3 is the modified skewness of the signal: 



AM 

A-0 



4"^ Older moment^m4rwcalculated.-m4 is the modified lairtosis-Of the 

20 signal: 



JV-l 
a«0 



2.5 Data analysis in Matlab 



25 The T-wave moiphology parameters for the acquired, pre-processed ECQ recordings 
were evaluated using Matlab 6.0. Only valid data were analyzed - i.e. data from leads 
where the signal vims not conupted by high frequency noise and where the evait detec- 
tion algorithm was successful in detecting the relevant events with satisgictory pred- 
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ston. Parametei means and srandaid deviations were calculated fcr every T-wave in 
die signal on all leads. A great interlead variation in T-mve morphology may be an 
indicator of WJTS. Interlead variance was therefore examined by calculating the stan- 
ch deviation of the lead means for each paiamelra. 

■ 

5 

Only the parameter means from lead V5 and interlead standard deviations were used 
as final parameter values. Hence, fbr every parameter in table 1, two paiameteis were 
calculated - one with index "meanVS" and one with index «std" e.g. FlmcanVS and 
Flsid. 

10 2.7 Statistical analyste 

In order to characterize and classic data from die three groups CKvLQTl, HERO and 
normal), the evaluated parameter values were processed using disciinunant analysis. 
Hie analysis was carried out in SPSS version U.5. The objective of the discriminant 
analysis was twofold: finding parameters that most efBcaently discriminate 

15 the groups and reducing the number of variables. Therefore astep«wise procedure was 
used with the Mahalanobis D2 as the most appropriaw distance measure. 

The enny/temoval-criteria vrere adjusted in order to reduce tiie nnmber of variables in 
the discriminant fimctions to achieve a 1 :3 ratio between the number of variables and 
20 the popuJatiOD size (N=l 6)- Ths criteria were empirically chosen to be p«itry - 0.04S 
and prcmoval = 0.09 providing the desired 5 variables in the discriminant fimctions. 

Figure 4, Scatterplot showing classification of individuals by genotype. Separation of 
groups was carried out by 2 discriminant fimdiflns vAib. 5 variables Aat diaiacierize 
25 r^olarisation by computation of symmetry, flatness and duration. 

3, ResuUs 

Tha discriminant functions were based on data from aU KvLQIl , HERO and normal 
subjects. TUB 5 parameters included in both discriminant functions arc listed in table 2. 
The discriminative efBdency of both generated functions was statistically significant 
30 after inclusion of aU 5 parameters (function 1: p<Q.0001. function 2: p<0-005). 
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Variables Entered 



Step 


Entered 


1 


Fllstd 


2 


QTcmeanVS 


3 


S5meanV5 


4 


D4std 


5 


S4ineanV5 



Table 2. Variables used by the two discriminaUng functions. Stepwise introduction of 
more variables improved Hie abUity of the functions to discriminaie between KvLQTl. 



HERO and normal. 



A scatteiplot was generated from the discrimination functions and groupings of indi- 
10 vidual genotypes can be seen in figure 4. The dotted lines were read fi-om the SPSS 
generated terriiorial map and manually added. The lines reflect borderlines where the 
differences between each pair of discrimination functions are zero. All 16 processed 
ECG's were correctly classified and showed at least one discriminatoiy characteristic 
as defined by the 5 parameters included in the discrimination functions. Cross valid*- 
15 tion of bodi discriminant functions was done with the leave-one-out method and aU 1 6 
subjects were again correctly grouped. Reducing the number of variables resulted in 
misclassified cases due to lack of one or more discriminatory characteristics. In light 
of this finding we elected to pexform further analysis of the selected parameters in or- 
der to investigate the individual contributions of each variable to the separation of the 
20 tlwee primary groups of subjects. Extreme values for all parameters were identified 
and the mean was computed. 

The result is plotted in figure 5. As expected the extent of interiead flatness variation 
observed in HERO and normal individuals was lower than that found ia KvLQTl sub- 
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jects. This is evidenced the Fl Istd parameter in figure 5a. When evaluating pa- 
lameter values S4meanV5 and S5meaaV5 (figure 5b» d) the extent of 6/8 asymmetry 
in KvLQTl and nonnal was generally less than tbai of HERG individuals. Both 
S4meanVS and SSroeanVS ai« synunetry paramoters and asynunetiy in HERG indi- 
viduals was augmented in two ways: When bifid T-waves were present the interval 
from Tsian to Ttop was prolonged due to the definition of Ttop used in this study (the 
last highest point on the T-wave). Also, when the initial portion before Ttop was pro- 
longed in HERG individuals better discrimination was possible. Both phenomena 
were observed in HERG subjects. Generally the Ba2!etr corrected QTc observed in 
HERG and KvLQTl was higher than that of normal mdividuals (figure 5e). However 
overiap existed between all three groups preventing separation of the groups by QTc. 
Since no single parameter included in the discrimination fimctioos was able to sepa- 
rate KvLQTl. HERG and normal, we proceeded to investigate the classification effi- 
ciency provided by the three primary categories represented by the parameters in the 
fitoctions. This was carried out by generating new discrimination fimctions using pa- 
. rameteis fiom one category only while e.KcIuding the other two. Then, from the new 
discrimination fiinctions three additional fimctions were generated, this tinje allowing 
the inclusion of parameters from corabinations of two categories. Scattetplots illustrat- 

* _ 

ing the results of this analysis are shown in figures 6a-f. The first two functions (figure 
6a) included parameters that characterize the symmetrical properties of the Twave. 
83.1% of the 16 subjects were correctly daasified. Arrows in figure 6a indicate the 3 
raisdassified subjects. A second discriminant analysis was perfonned using flamess 
parameters. This resulted in 93.8% correctly classified subjects. Only one subject was 
not correcUy classified as indicated by the anow on figure 6b, The misclassified case 
was the same HERG subject inconecfly classified ushag symmetry parameters. The 
discriminatory efficiency of duration parameters was also evaluated. Discrimination 
analysis resulted in 93.8% conrectly classified subjects. One HERG subject was mis- 
classified as KvLQTl. QTc was 416ms and the subject showed relatively peaked T- 
waves similar to those found in KvLQTl. However the duration parameters fiiiled to 
identify this morphological feature, thus reducing classification performance. 

It can be noted that improved classification was obtained using flatness or duration 
parameters versus symmetry parameters and it seemed reasonable to investigate if fiir- 
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ther classification improvement could be adueved using a combination of several pa*- 
rameeer categories. 

Figure 5, a) Fllstd -Interlead standard deviation of t&e ratio between the total area 

* 

5 under tlie T-wavc fyom Ttop to Tend and the coiresponding time interval, b) 
S4meanV5 - Lead VS mean modified skewness evaluated in a symmetrical interval 
surrounding Ttop and corresponding to 20% of the interval between Tstait-Tend. c) 
D4$td - Interlead standard deviation of the time interval from Ttop to Tend, d) 
SSmeanVS - Lead V5 mean of the ratio between the time interval from Tstart to Ttop 
10 - and the corresponding time interval from Ttop to Tend, e) Lead V5 mean QTc, 

Figures 6d-f show the results of three separate discritninant analysis using combina- 
tions of parameters from two categories. It can be noted that classification of subjects 
was perfect in all cases, even when repolarisation duration was not considered (figure 
15 6d). 

4. ConclusioA and discussion 

J* 

Tlic initia] discriminant analysis performed in this study resulted in perfect classifica- 
tion of all KvLQTl, HERO and nomial subjects. In table 2 it was norcd that ihe dis- 
criminant functions included parameters from all three categories; T-wave symmetry, 
20 T-wave flatness and duration. This is in agreement with the initial hypothesis that a 
combination of repolarisation duration and T-wave morphology characteristics could 
improve discrimination between KvLQTl, HERG and normaL 

To xmderstand v/tiy some subjects were misclassified using a reduced set of parameter 
categories (figures 6a-c) the duration parameters and morphological characteristics of 
25 all 16 ECQ's were examined. 

Using only symmetry parameters, 3 subjects were misclassified. However no obvious 
visual characteristics on the three misclassified ECCs could be identified that ex- 
plained the incorrect classifications. The Bazetl corrected QTc was 347ms for the 
normal subject, 425ms KvLQTl « 476ms HERG. Although an obvioiisly prolonged 
30 QTc was present in the misclassified HERO subject it was not identified using sym- 
metiy parameters alone. 
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Discriminant analysis using parameters from the flatness category resulted in only 1 
misclassification. Again no visual characteristics were identified to account for the 
misclassificaiion. Although it was anticipated that the 

Figure 6, a) Tbe result of discriminant analysis using symmetry parameters resulted in 
. Jee misclassified cases (arrows). Visual inspection of the ECG's revealed no appar- 
ent abnormalities to indicate the reason for incorrect misclassification. b) Hie result of 
discriminant analysis using flatness parameters. One inootrectly classified HERG sub- 
ject was identified (anow) even though no obvious visual abnormality indicated a dif- 
ferent genotype, e) Result of discriminant analysis using duration parameters. This 
result iUustraies the feilute of duration parameters to discriminate between KvLQTI, 
HERG and normal (arrow), d-e) Combinations of paiaraetets from two categories il- 
lustrate the improvement in classification efficiency when compared to figures 6a-c 
evaluation of T-wave flamess would be able to discriminate HERG from KvLQTl 
subjects this was not accompUshed by using flatness as a single descriptor of repolari- 
sation. Performing discriminant analysis based on the QTc parameter as the only vari- 
able tesuhed in 1 misclassification- This was not unexpected smoe it is weU known 
that a substantial overlap in QTc values can eydst between noimal and affected indi- 
viduds. The lack of unambiguous discrimination between all groups by use of the 
QTc parameter alone emphasizes the hypothesis diat additional parameters are needed 
to classify LQTS individuals. By combining parameters fiom two categories it was 
found that the discriminatory strength was hiweased. 



(figures 6d-0 This was evidenced by the fact that no suDjccis were huscbu.^.^^. 
two categories. A pardculariy interesting finding, was the perfect separation of all sub- 
jects that was obtained using symmetry and flamess parameters with no duration pa- 
rameters included. This result implies the discrimhiatory strength inherent in parame- 
ters fiom those two categories. In addrtion it was fmjnd that symmetry or flatness pa. 
lamfeters combined with duration parameters yielded perfect discrimination between 
all groups. Results from the discriminant analysis using One and two categories indi- 
cate that a combination of more parameter categories strengthen the overall discrimi- 
natory power of the classification fimctions. Combining these findings with the results 
from the three category discriminant analysis initially performed, it is reasonable to 
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speculate that a substantially improved discrimination between KvLQTl, HERG and 
nonnal is possible using all tlwee categories of parameters. 

In light of the results obtained in this study we propose a new technique for discrimi- 
5 nating between KvLQTl, HERG and normal subjects. Through multivariate discrimi- 
nant analysis it was found that a combination of two duration parameters and three T- 

« 

wave symmetiy-and flatness parameters was sufficient to classify each of the 16 study 
subjects into one of the three distinct groups. Although no single parameter had the 
necessaiy discriminatory strength to classify the subjects, die combination of multiple 
10 parameters in two discrimination fimctions was statistically significant (function 1; 
p<0.0001, function 2: p<0.005). The encouraging results of multivariatie repolaiisation 
analysis found in this study support the use of syrametiy-, flatness- and duration pa- 
rameters to classify LQTS patients. 

IS Furflier effort must be made to strengthen the statistical impact of this study and to 
investigate whether the proposed synunetry-, flatness- and duration parameters are the 
best suitable for the discriminatioa functions. 



The use of the proposed multiple parameter categories to classify KvLQTl and HERO 
20 genotypes may prove to be a powerful clinical tool in the making. 



25 



The invention will in the following be described in detail with reference to the draw- 
ing, where figures 1 to 6 are mentioned in the text, and they are as such not furth^ 
described. 



Fig. 7 shows examples of £CG-curvature having ST-elevation myocardial infarction, 
where the active parameters for indicating ST-elevation myocardial infarction is ST- 
elevation, ST-morphology, T-wave morphology and Q-wave morphology. 

30 Fig, S shows curvaturic non ST-elevation myocardial infarction, where parameters of 
ST depression, T-wave morphology or Q-wave moiphology could be used. 
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Fig, 9 shows a ECG-curvature indicating Cardiomyopathia where the following pa- 
mmeters can be used for indication: P-wavc morphology, QRS duration, S-wave mor* 
phology and T-'wave tnoiphology. 

5 

Fig. 10 shows in the same way indications for Biugada Syndrome. Here, tive effective 
paiameters for indicatioa are PR-duration, ST-elevation, ST motpholo^ and T-wave 
moiphology. 

■ 

10 Fig. 1 1 shows curvatures referring to Right bundle branch block REBB vtdiere the pa- 
rameters could be QRS duration, R-wave moiphology, T-wave morphology and ST- 
elevation. 

Fig. 12 shows curvatures indicating Left bundle branch block LBBB where the effec- 
1 5 tive parameters could be QRS duration, R-wave morphology and T-wave morphology. 

Fig, 13 shows curvatures indicating Short QT syndrome where parameters for indica- 
tion could be Q-T duration and T-wave morphology. 

* 

20 Fig. 14 shows curvatures indicating Hyperkalemia where parameters efifective for in- 
dications are P-wave moiphology, T-wave morphology^ QRS duration^ QT duration 

■ 

and PR duration. 

Fig. 15 shows curvatures indicating Hypokalemia where the effective parameters for 
25 indication seem to be QT duration, T-wave morphology and ST depression. 

Fig 16 shows curvatures indicating Pericarditis where the effective parameters are ST 
elevation, ST morphology, Q-wave morphology and PR depression. 

30 Fig. 1 7 shows a curvature indicating Right Ventricular Hypertrophy (RVH) where tlie 
effective parameters for indication are Q-wave morphology, QRS duration, S-wave 

morphology and T«wave morphology. 

« 

* 
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Fig. 18 shows curvawes indicating Left Ventricular Hypertrophy (LVH) where the 
effective parameters for indication are Q-wave morphology. QRS duration, S-xvave 
morphology and T-wave morphology. 



5 Fig. 19 shows curvatures indicating Arrhythmogenic Right Ventricular Dysplasia 
where the parameters are QRS duration, S-wave moiphology and T-wave morphology 
is used. 
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CLAIMS 

1 . A system for analysing ECG cxavatute vdierein at least one among a number of 
diffezoxt parameters is isolated, which system has input means connected to an ECG 
source, where the different parameters of a received ECG curwrture are indicated 
and/or isolated for indicating possible symptoms which relates to or are indications of 
certain deceases, where said deceases are knovai to influence the ECG curvature, and 
that a first number of selected parameters, are combined in at least a first mathematical 
analyse, where tlie result of flie analysis is represented as a point in a cocxdinate sys- 
tem, comprising at least two axes, where the system compares the actual placement in 
the coordinate syst^ with a number of reference parameters stored in the system, for 
indicating symptoms of diseases having influence on the ECG curvature, 

2. System according to claim I,charactcri aedinthattheanalysingpnocess 
is repeated in the system for further selected parameters in order to for indicating fur- 

m 

' ther possible diseases or symptoms. 

3. System according to claim 1 or 2, character fsedin that the parameters 
are divided into at least two main groups in the system, which groups contains pa- 
rameters of synunetiy, flatness and duration. 

4. System according to claim 3. characterized in that the group of symmetry contains 
at least following parameters: 

5 1 Skewness evaluated from Tstart to Tend, 

■ 

52 Skewness evaluated from Tsiart-to Tend with Ttop as mean, 

53 Skewness evaluated in a symmetric interval, 10 % of the Tstart-Tend interval 
surrounding Ttop widi Ttop as mean 

54 Skewness evaluated i a symnietric interval, 20 % of the Tstart-Tcnd inierval 
suTTOunding Ttop wit Ttop as mean, 

85 Ratio of the time interval from Tstart to Ttop and the time interval from Ttop to 
Tend. 
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S6 Ratio of the average slope from Tstart to Ttcp and from Ttop to Tend, 
where the group of flatness contains at least the following parameters: 

F 1 KurtosLs evaluated from Tslart to Tend, 

F2 Fl nonnalized by the absolute Rtop-Qnadir value, 

F3 Kurtosis evaluated from Tstart to Tend with Ttop as mean, 

F4 F3 nonnalized by absolute Rtop-Qnadir vaJue^ 

F5 Kurtosis evaluated in a symmetric interval, 10 % of tlie Tstart-Tend interval sur- 
rounding Ttop with Ttop as mean, 

F6 F5 normalised by absolute Rtop-Qnadir value, 

F7 Kurtosis evaluated in a syxnmetric interval, 20 % of the Tstart-Tend interval sur- 
rounding Trop with Ttop as mean» 

F8 Kurtosis nonnalized by the value of Rtop with Ttop as mean, 

F9 Ratio of the total area under the T-wave from Tstart to Ttop and the correspond- 
ing time interval, 

FIO. F9 nonnalized by absolute Rtop-Qnadir value, 

Fl 1 Ratio of the total area under the T-wave from Ttop to Tend and the correspond* 

/ 

ing time interval, 
F12 Fl 1 normalized by absolute Rtop-Qnadir value 

F13 Ratio of the total area under the T-wave from Tstart to Tend and the correspKDnd* 
ing time interval, 

F 14 F 1 3 normalized by absolute Rtop-Qnadir value, 

F15 Ratio of the height of Rtop and the width of the Tstart-Tend interval, 

where the group of duration contains at least the follovsdng parameters: 

QTc The Q*T interval nonnalized by the square root of the R-R interval according to 
Ba^ett^s formula, 

D2 . Time interval from Tstart to Tend, 
D3 Time interval from Tstart to Ttop, 
D4 Time interval from Ttop to Tend. 
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5. System according to claim 3, c h a r a e t e r I s c d in that the system arc com- 
bining panamet^ from di£fer^t groups. 

6. System according to one of the claims l-4,characterSscdmthatthesy&- 
lem analyses the QT interval of the ECG curvature for indicating Long QT syndrome. 

7. Method for analysing ECG curvature, which curvature contains a number of pa^ 
rameters, c h a r a c t e r i s c d In that the method for analysing tiie ECG curvature 
incorporates the steps of: 

a) receiving ECG curvatuie from a source, 

b> indicaie a number of different parameters contained in the received ECG curva- 
ture, 

c) storing the parameters in storage means, 

d) select disease specific parameters in the storage means 

e) combine selected parameters in mathematical analysing means 

f) represent the result of the mathematical analyse as a point in a coordinate sys- 
tem, comprising at least two axes» 

g) compare the actual placement in the coordinate system with a number of refer- 
ence parameters stored in a memory, 

h) indicating diseases having influence on the ECG curvature. 

8, Method according to claim l8^-cJi-a r a c t e r i s e d in that the method is re- 
peating the analysing process ftit Jastiha selected parameters for indicating fxjrther 
diseases. 

9. Method according to claim 18 or.l9, c h a r a c t c r i s e d in that the method are 
dividing parameters into at least three groups, which groups contains parameters of 
symmetry, flatness and duration. 



* 
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ABSTRACT 

The present invention relates to a system or a method for analysing ECO curvature 
where at least one among a number of different parameters is isolated^ which system 
has a input means connected to an ECG source, v^exe the different pacameters of a 
5 received ECG curvature are indicated and/or isolated and for indicating possible 
symptoms which relate to or are indications of certain deceases, where said deceases 
are known to influence the ECG curvature. 

The aim of the invention is to achieve a system and a method for objective, fast and 
10 effective indication of a number of symptoms derivable from an ECG curve which 
may be indicative of one or more diseases. 

This can be achieved with the system or method previously described, if a first number 
of selected parameters are combined in at least a first mathemadcal analysis, where the 

IS result of the analysis can be represented as a point in a coordinate system comprising 
at least two axes where the system can compare the actual placement in the coordinate 
system with a number of reference parameters stored in the syst^ for indicating dis« 
eases having influence on the ECG curvature. Hereby, it is achieved that any symptom 
of a disease having an indication (influence) in the ECG curvature can be detected in 

20 an objective, automated and very fast way. The system might be used under field con- 
ditions such as in ambulances or in other situations where a fast indication of heart 
diseases is needed in order to help the patient in a correct way as early as possible. 

* 

Fig. 4 
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Function 1 
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Fig. 5b 
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Fig. 5c 
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Fig. 5d 
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Fig. 5e 
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Fig. 6b 
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Fig. 6c 
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Fig. 6d 
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Fig. 6e 
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Fig. 6f 
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ST-elevation myocardial infarction 




• ST morphology 

. • T-wave morphology 

• Q-wave morphology 

Traditionally the presence of an ST elevation myocardial infarctioxi (STEM!) is in 
iht ECQ diagnosed in tbe early phase by the ST elevation In the J-point and the 
subjective evaluation Of the shape of the ST-moxphology cogeflier with tbe 
progressive development of Q-waves. With the present method the diagnosis wfll 
be based on tbe combination of d)e objective measures and new parametcis inakiiig 
the subjective changes into objective potameters. E.g the shape of the ST 
segment (ST movphology) will be characterised by 1 or more parameters, by which 
it is expected that it wiU be possible to iinprove the diagnosis in an automated way. 
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Fig. 8 



Non ST-elevation myocardial infarctioi 




• • ST depression 

• T-wave morphology 

• Q-wave morphology 
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Fig. 9 



Cardiomyopatliia 
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P-wave morphology 
QRS duration 
S-wavc morphotogy 

T-wave morphology 
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Fig. 10 



Brugada Syndrome 
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• PR-duration 

• ST-elevation 
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• ST morphotogy 

• T-wave morphology 



EfliPfaneszd j t IS-Dez* 1 1 :30 



Fig. 1 1 




• QRS durntiaii 

• R-wave morphology 
» T-wave morphology 

• ST-elevation 



Fig. 1 2 

Left bundle branch block LBBB 




t QRS duration 

• R-wave morphology 

• T-wave morphology 
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Short QT Syndrome 



• mm m. « J ^ ^ 




' Mr 



!^':Biff^ 




Kb 



B 



Q-T duration 
T-wave morphology 



EiDPfansszei t Id-Dez* 



11:30 



Fig. 14 



Hyperkalemia 




A 



— ./\. 



Chan^fi^ In iSie ECG in taad 0 
caused by hyipefkalamla 



• P-wave morphology 

• T-wave morphology 

• QR5 duration 
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Fig. 15 



Hypokalemia 
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Changes tn the ECG irt tead It 
caused by hypokalemia 
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Fig. 16 



Pericarditis 
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Bwlfln «rty tvppteiisilon ^ ST/ T Acute pctfaxtRCis wUh STI T 



• ST elevation 

• ST morphology 

• Q-ivave morpbology 

• PR depression 
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Right Ventricular Hypertrophy (RVH) 
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Left Ventricular Hypertrophy (LVH) 
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Q-wave morphology 
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Fig, 19 



Arrhythmogenic Rigkt Ventricular Dysplasia 




• QRS duration 

« S-wave morphology 

• T-wave morphology 



EeiPf ansszei t 19.Dez. II :30 



PCT/DK2004/000722 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record. 



Defective images within this document are accurate representations of the original 
documents submitted by the apphcant. 

Defects in the images include but are not limited to the items checked: 

^^M.ACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 



□ OTHER: ^ 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



BEST AVAILABLE IMAGES 




LINES OR MARKS ON ORIGINAL DOCUMENT 



REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 



